General analytic formulae for attractor solutions of scalar-field dark energy models 

and their multi-field generalizations 
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We study general properties of attractors for scalar-field dark energy scenarios which possess 
cosmological scaling solutions. In all such models there exists a scalar-field dominant solution with 
an energy fraction = 1 together with a scaling solution. A general analytic formula is given to 
derive fixed points relevant to dark energy coupled to dark matter. We investigate the stability of 
fixed points without specifying the models of dark energy in the presence of non-relativistic dark 
matter and provide a general proof that a non-phantom scalar-field dominant solution is unstable 
when a stable scaling solution exists in the region < 1. A phantom scalar-field dominant fixed 
point is found to be classically stable. We also generalize the analysis to the case of multiple scalar 
fields and show that for a non-phantom scalar field assisted acceleration always occurs for all scalar- 
field models which have scaling solutions. For a phantom field the equation of state approaches that 
of cosmological constant as we add more scalar fields. 

PACS numbers: 98.80.-k 
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I. INTRODUCTION 

Dark energy is one of the most serious mysteries in 
modern cosmology. The observations of supernova (SN) 
la have continuously confirmed that about 70% of 
the total energy density of universe consists of unknown 
energy which leads to an accelerated expansion of the 
universe. The existence of dark energy has been also sup- 
ported from the observations of cosmic microwave back- 
ground (CMB) Q. large-scale structure (LSS) and 
baryon oscillation experiments (see also Refs. 

While cosmological constant can be a natural can- 
didate of dark energy, this suffers from a severe fine- 
tuningproblem if it originates from particle physics (see 
Refs. IE S B 13 for review). Instead of sticking to cos- 
mological constant, a host of alternative scenarios have 
been proposed- ranging from scalar-field models, Chap- 
lygin gas [T^j, braneworld Hi) and modified gravity sce- 
narios in General Relativity [TJ- Since scalar fields are 
natural ingredients in particle physics, it is an interest- 
ing challenge for theoreticians to construct viable dark 
energy models using them. A number of scalar-field mod- 
els have been studied as candidates of dark energy- an 
incomplete list includes quintessence , k-essence [l4| , 
tachyons E^l, phantoms [Tt| and (dilatonic) ghost con- 
densates [13, [3 • These can at least alleviate the fine- 
tuning problem of cosmological constant because of their 
dynamical nature. 

A general guidance to obtain viable dark energy mod- 
els is that the energy density of it is "hidden" during 
radiation and matter dominant eras and comes out near 
to the present epoch. If one uses the pro perty of cos- 
mological scaling solutions 0, |2jJ, l2ll l22j. the energy 
density of scalar fields decreases proportionally to that 
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of the background fluid (radiation or matter) indepen- 
dently of the initial conditions [2j| . The system can exit 
from the scaling regime to give rise to the present ac- 
celerated expansion if the potential of the field becomes 
shallow [23, [24| °r if the coupling Q between dark energy 
and dark matter becomes important [2f|. Especially in 
the latter case there is an intriguing possibility that the 
present universe is a scaling attractor with a constant 
dark energy fraction £1$ ~ 0.7 [2(| ■ Thus the presence of 
scaling solutions plays an important role for the model- 
building of dark energy. 

It is well known that exponential potentials give rise 
to scaling solutions for an ordinary scalar field |2J|. The 
form of the potentials which generate scaling solutions 
depends on the models of dark energy. For example in 
the case of a tachyon field the corresponding potential is 
given by V{4>) = Vo^ 2 [H |H |29| . A S eneral algorithm 
to obtain potentials corresponding to scaling solutions 
is given in Refs. [Tsl . l30l |. Interestingly the scalar-field 
Lagrangian density for the existence of scaling solutions 
is restricted to be in a compact form p = Xg(Xe x ^) [l8j . 
where X is the kinematic term of a scalar field <f> and g is 
an arbitrary function. This analysis can be extended to 
the more general case in which the Hubble parameter H 
has a dependence H 2 cx p n with p being the total energy 
density |30j. By using the formula in Ref. |30j we can 
easily obtain potentials which generate scaling solutions 
once we specify dark energy models. 

In addition to scaling solutions there exists another 
important fixed point which corresponds to the scalar- 
field dominant solution characterized by 0,$ = 1. In fact, 
in the case of an ordinary scalar field with an exponential 
potential, this solution was used in Ref. [23| to obtain a 
late-time attractor with an accelerated expansion after 
the system exits from a scaling re gim e during radiation 
and matter dominant eras. In Ref. ml it was shown that 
the scalar-field dominant solution also exists for tachyons 
and dilatonic ghost condensates. In this paper we will 



2 



show that this fixed point always exists for all scalar- 
field models which possess scaling solutions. We shall 
also provide a general analytic formula to derive critical 
points for both the scalar-field dominant solution and the 
scaling solution. 

The stability of critical points against perturbations 
was analyzed in Ref. pul ] by specifying three models of 
dark energy. We shall show that the stability analysis 
can be analytically done for all dark energy models which 
have scaling solutions. We will generally prove that the 
non-phantom scalar-field dominant solution is unstable 
when the stable scaling solution exists in the region < 
1. It will be also shown that the phantom scalar-field 
dominant solution is classically stable. This analysis is 
useful when we construct viable dark energy models. 

There is another interesting possibility to give rise to 
an inflationary solution by using multiple scalar fields. 
In the case of an ordinary scalar field with an exponen- 
tial potential many fields can cooperate to sustain infla- 
tion even if none is able to do individually This 
assisted inflation scenario was applied to dark energy in 
Refs. [33|,li3- The field energy density mimics that of the 
background fluid during radiation and matter dominant 
eras because of the scaling property, which is followed by 
the scalar-field dominant solution as more fields join to 
sustain an accelerated expansion. In this paper we shall 
show that for a non-phantom field this assisted behavior 
always happens for all scalar-field dark energy models 
which have scaling solutions. 



II. CONDITIONS FOR THE EXISTENCE OF 
SCALING SOLUTIONS 

In this section we shall review the Lagrangian for the 
existence of scaling solutions presented in Refs. 0, |3(J • 
Let us start with the following 4-dimcnsional action 



S 



R 



(i) 



where R is a scalar curvature and X = ~g fJ, ' / d fl (j)d v (j>/2 is 
a kinematic term of a scalar field (f>. S m is an action for a 
barotropic fluid which is generally dependent on <f>. The 
reduced Planck mass M p \ is set to be unity. The above 
action covers most of scalar-field dark energy models. 

The pressure density and the energy density of the field 
are given, respectively, by 



P<p=P, P<t> = 2X ~Qx~ p - 
We introduce a scalar charge a defined by 

1 5S m {^) 



(2) 



(3) 



which corresponds to the coupling between the field <ft and 
the barotropic fluid. This coupling can lead to an accel- 
erated expansion even when scaling solutions can not do 



so in the absence of the interaction |25j . See Refs. |35j 
for recent works about the interacting dark en ergy . A 
similar coupling arises in scalar-tensor theories [25l Isf3 | 
and in neutrino models of dark energy in which the mass 
of neutrinos depends on a scalar field <j) EI3 • 

In a flat Fricdmann-Robertson- Walker (FRW) back- 
ground with a scale factor a, the continuity equation for 
the field <fi is given by 



3H(1 + w^p^ = -Qp r , 



(4) 



where a dot represents a derivative with respect to cosmic 
time t, Q = u/p m , = p$/p<$> and H = a/a is the 
Hubble rate. Meanwhile the fluid energy density, p m , 
satisfies 



p TO + 3H(l + w m )p m = Qp m (j> . 



(5) 



In what follows the equation of state of a barotropic fluid, 
w m = p m / p m , is assumed to be constant. 

The Hubble parameter H obeys the Friedmann equa- 
tion 

3H 2 = p<j> + p m . (6) 
The fractional densities of p§ and p m are defined by 

(7) 



o = ^ o - Pm 



These satisfy the relation + Q m = 1 from Eq. © . We 
note that reconstruction equations suitable for the com- 
parison to SN la observations were obtained in Ref. [3^ 
for the action 

Cosmological scaling solutions are characterized by the 
condition 



P<f,/pm = const 0) 



(8) 



In this case fl^ and O m are constant from Eq. (JZJ). Wc 
also assume that Q and do not vary in the scaling 
regime. Note that the constant coupling Q appears in 
theories in which the mass of dark matter or neutri- 
nos depends exponentially on the scalar field, see e.g., 
Refs. [Hill. 

From the condition JHJ wc find (logp^)* = (logp m )V 
Then we obtain the following relation from Eqs. and 
©: 



dTV 



(w m -w^) = const . 



(9) 



where N = J Hdt is the number of e-foldings. This gives 
the following scaling behavior: 



dlnp^ dlnp,- 



dTV dTV 
where the effective equation of state is 
w^ptp + w m p m 



3(1 + W e ft) , 



(10) 



Wcff 



P4> + Pm 



W m + Q^W^ ~ W m ) . (11) 



3 



This expression of w e g can be used irrespective of the 
fact that scaling solutions exist or not. 
The definition of X gives the relation 



2X = H ( d^J K H <X(P* + Pm), (12) 



field with an exponential potential: p = X — ce x ^ 
HI 113. The function g(Y) = -1 + cY gives the dila- 
tonic ghost condensate model with Lagrangian density 
p = —X + cX 2 When we choose the function 

g(Y) = — cVl — 2Y" jY together with the introduction 
of a new field ip = (2/A)e A(?i / 2 , we obtain the tachyon 



Lagrangian p = —V((p)y/T~- 
where we used Eq. Hence the scaling property of X potential V((p) — Ac/(X 2 ip 2 ) 
is the same as and p m , which leads to the relation 



- Lp 2 with an inverse square 
I2II2I. 



dlnX 
~~dN 



= -3(1 + W e s) . 



(13) 



The pressure density p^ = w^p^ scales in the same way 
as p0, i.e., dlnp^/dTV = — 3(l + w e ff)- Using this relation 
together with Eqs. © and lfT3)l we find 



Qlnp^ l^lnp^ 



d\nX X 



1 



where 



A = Q 



1 + w. n 



Q,<h(w ri 



w 4>) 



(14) 



(15) 



From Eq. (|14|> the existence of scaling solutions re- 
stricts the Lagrangian density in the form: 



p{X,<l>)=Xg{Xe x *) , 



(16) 



where g is an arbitrary function in terms of Y = Xe^. 
This expression was first derived in Ref. |l8j and was 
extended to the more general case in which the Fricd- 
mann equation is given by H 2 oc (p^ + p m ) n |3Cj. The 
quantity Y is conserved along the scaling solution, i.e., 
Y = Xe Xct > = constant. 

Although we derived the expression (|16f) under the as- 
sumption that Q is a constant, this can be generalized 
to the case in which Q depends on the field <f>. This case 
will be presented in elsewhere ^3 ■ We also note that the 
r.h.s. of Eq. © is kept to be a constant for Q — ► 0, since 
Wfj, — > w m for scaling solutions in this limit. Thus the 
above discussion for the derivation of Eq. I|16|) is valid 
even for Q = by taking the limit Q — > 0. 

Equations @, i|12[l and (|15|l give the relation 



3H 2 



2(Q + A) 2 



3(1 + w m ) 2 
From Eqs. (fTTf t and ltT5|) we find 



X. 



A-Q 



Q + A 



(17) 



(18) 



When Q = this reduces to w e s = w m . Hence an ac- 
celerated expansion is not possible for a non-relativistic 
fluid (w m = 0). However the presence of the coupling 
Q leads to an accelerated expansion (iu e g < —1/3) for 
Q > (3w m + l)A/2. 

If we choose the function g(Y) = 1 — c/Y in Eq. (|16fl . 
we obtain the Lagrangian density for an ordinary scalar 



III. FIXED POINTS 

In the previous section we showed that the Lagrangian 
(|16fl possesses scaling solutions. There exist other fixed 
points for the system (|16fl . We shall derive fixed points 
relevant to dark energy and analyze the stability of them 
against perturbations. Note that critical points were de- 
rived for three dark energy models in Ref. [3lT | by speci- 
fying the functional form of g(Y). In this paper we show 
that this is possible for all scalar-field models that pos- 
sess scaling solutions. In what follows we shall consider 
the case of positive values of Q and A. 

From Eqs. (0) and |@} together with the Lagrangian 
density (|16fl . we find that the scalar field obeys the equa- 
tion of motion: 

+ 'iHAp x 4> + XX [1 - A(g + 2Yg')] = -AQp m . (19) 

where a prime represents a derivative with respect to Y 
and 

A(Y) ee [g(Y) + 5Yg'(Y) + 2Y 2 g" {Y)]^ , (20) 
Px(Y)^^=g(Y)+Ygi(Y). (21) 
Combining Eqs. (0J and JSJ with (JSJ, we find 

H = -X Px -^(l + w m ) Pm . (22) 

We introduce the following dimcnsionless quantities 
which are useful to study the dynamical system: 



-A0/2 



V6H 



V3H 



(23) 



The variables x and y can be regarded as a "kinematic" 
term and a "potential" term, respectively. Since we do 
not consider a contracting universe, y is positive from its 
definition. The variable Y = Ae^ is expressed in terms 
of x and y: 



Y = x 2 /y 2 



(24) 



Equation © gives the constraint equation 0^+f2 TO = 1 
with 

n 4> = x 2 (g + 2Yg'). (25) 
The equation of state for the scalar field <fi is given by 

9 



W<f, 



9 + 2Yg' 



(26) 
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It may be useful to notice the relation 



Scalar- field dominant solutions 



2x 2 

= gx 2 , w t f > =—l + — — px ■ (27) 
11,/) 

This shows that the field behaves as a phantom (w<f, < 
— 1) for px < provided that VL^ > 0. 

From Eqs. (|19|l and (|22|l we get the following au- 
tonomous equations for x and y: 



d.r 
ON 



dy_ 

dN 



3x 



1 + gx 2 - w m {$l^ - 1) 
V6A 



X.r 



(Q + X)n^-Q-V6(g + Yg') 



3y 
2 



1 + gx 2 - w m (^0 - 1) 



V6 



Ax 



.(28) 
(29) 



We can obtain fixed points of the system by setting 
dx/dN = and dy/dN = 0. From Eq. © we find 
that y = corresponds to one of the fixed points. How- 
ever this is irrelevant to dark energy, since we need the 
contribution of a "potential" term to give rise to an accel- 
erated expansion. In addition this fixed point is unstable 
against perturbations |31j . 

Then the fixed points we are interested in satisfy the 
following equations: 

V6Xx = 3 [l + gx 2 - w m (ft - 1)] , (30) 
V6{g + Yg')x = (Q + X)fl^ - Q . (31) 

Using Eq. ^Jj we find that g + Yg' = f^(l + w <t> )/2x 2 . 
Then Eqs. I|3U|) and l|31|l can be written in the form: 



V6[l + (w<j, - W m )Q<f, + W n 

2A 

VG(1 + w^Q^ 



(32) 
(33) 



2[(Q + A)Q -Q] ' 
which leads to 
(fy - 1) [(w+ - w m ){Q + A)ft + Q(l + w m )} = . (34) 
Hence we find two cases: 
• (i) A scalar-field dominant solution with 

= i . (35) 



• (ii) A scaling solution with 
(1 + w m )Q 



fid 



(w m - w<f,)(Q + X) 



(36) 



In the case of an ordinary scalar field, for example, the 
solution (i) corresponds to the 4-th fixed point in Table 
I of Ref. [20|, whereas the solution (ii) corresponds to 
the 5-th fixed point. The above discussion shows that 
similar fixed points exist for all scalar-field dark models 
which possess scaling solutions. In what follows we shall 
study the properties of the critical points in more details. 



When 0^ = 1, Eqs. $27$ and (|3Hjl give the equation of 
state: 



V6A 

wj, = — H — —x . 



(37) 



In this case the effective equation of state of the system 
is given by w e s = by Eq. lllfl . The condition for an 
accelerated expansion, w c s < —1/3, corresponds to 



(38) 



From Eq. (|31J) we find that x is given by x = X/^/6px- 
Then uu is written in the form 



w<f> = -1 



3 P 



x 



(39) 



Hence one has w$ > — 1 for px > and w$ < — 1 
for px < 0. For an ordinary scalar field with an ex- 
ponential potential, the Lagrangian density is given by 
p = X — ce~ A< ^. Since px — 1 in this case the equation 
of state is w$ = — 1 + A 2 /3. This agrees with the result 
obtained in Ref. [2(| ■ In the case of a phantom field with 
an exponential potential (p = —X — ce~ A< ^), we obtain 
w <t> = -l- A 2 /3. 

From Eqs. lj3U|) and (j3T)l we find 



v6Ax — 3 
3^2 



Yg'(Y) = 



6- \Z6Aa 
6x 2 



(40) 



Once we specify a scalar-field dark energy model, i.e., the 
function g(Y), we get Y and x from Eq. I|40|l . Using the 
relation (|24J) we obtain the scalar-field dominant fixed 
point (x,y). The equation of state is known by using 
Eq. 123). 

Let us consider the ordinary (phantom) scalar field 
with g(Y) given by g(Y) = e — c/Y (e = +1 for the ordi- 
nary field and e = — 1 for the phantom). From Eq. H4()|l 
we find Y = cA 2 /(6 - eA 2 ) and 




(41) 



A 

V6e 



This agrees with what was obtained in Refs. |20tl31l| . 

In the case of dilatonic ghost condensate with g(Y) = 
-1 + cY, Eq. HI leads to 



Xl,2 

where 



V6A/ ± (A) 1 | A 2 / T (A) (42) 

4 ' cYl < 2 ~2 + ^^' (42) 



/±(A) = 1± N /1 + 16/(3A 2 ) 



(43) 



This result agrees with what was obtained in Ref. |18ll31j . 
We have two fixed points in this model. Since /+(A) > 



5 



and /-(A) < 0, the critical point (xi,yi) corresponds to 
a phantom one with < — 1 from Eq. (|37|) . whereas the 
point (£2)2/2) satisfies > — 1. In the latter case the 
condition H38[) for an accelerated expansion corresponds 
to A < \/6/3. The phantom divide is characterized by 
the condition p x = 0, i.e., cY = 1/2. From Eq. |g2J| 
we find 1/3 < cYi < 1/2 and cY 2 > 1/2. In the limit 
A — > both fixed points approach the phantom divide at 
C y = 1/2 (i.e., = -1). 



B. Scaling solutions 

When Qff, is given by l|5o)l. Eq. (J221 leads to 
_ _ \/6(l + w m ) 



2(Q + A) 



(44) 



This is equivalent to Eq. (|T7J> . From Eqs. ijTT]) and 

we obtain the effective equation of state given in 1)18(1. 

Equations (jjHJ an d (Bill give 



Q(<3 + A) + 3(1 + ui m )p 



A" 



(Q + A) 2 



(45) 



From Eqs. I|44|) and l|18fl we find that x and w; e ff are 
independent of the form of g(Y). Meanwhile Y and y 
depend upon the models of dark energy. From Eqs. i|30|) 
and we find 



(l-w m )g{Y)-2w m Yg'(Y) 



2Q(Q + X) 



(46) 



3(1 + w m ) 

Once the function g(Y) is specified, one gets Y and y = 



x/yY by using this equation. Then we obtain 0.$ and 
W,/, from Eqs. JJ^J and (|2"7ll . 

For example, in the case of an ordinary scalar field with 
g(Y) = 1 — c/Y, we obtain 



Y 



V 



3c(l + w m ) 2 



2Q(Q + A)+3(l-u&) 
2Q(Q + A) + 3(1 - u4; 



1/2 



(47) 
(48) 



2c(g + A) 2 

which agrees with the result obtained in Ref . [3l| . 

IV. STABILITIES OF FIXED POINTS 

In this section we shall study the stability of critical 
points obtained in the previous section. Remarkably this 
can be carried out without specifying the functional form 
of g(Y). Let us consider small perturbations 8x and 8y 
about the critical point (x Cl y c ), i.e., 

x = x c + Sx, y = y c + 5y, Y = Y c + 5Y . (49) 

At linear level 5Y is related to 5x and Sy via 



5Y = 2 l^Sx-^Sy 



(50) 



The function g(Y) is expanded as 



9(Y) = 9c + g' c (Y - Y c ) + ^(Y - Y c f + ■ ■ ■ , (51) 

where g c = g(Y c ). When we consider linear perturba- 
tions of g(Y), it is sufficient to neglect the terms higher 
than the first order. However when we evaluate linear 
perturbations of the g + Yg' term in Eq. I|28|) , we need to 
take into account the second-order term in Eq. Q51[l. i.e., 



S(g + Yg') = (2g' c + Y c g" c ) 8Y . 



(52) 



Similarly the perturbation of the fractional energy den- 
sity, = x 2 (g + 2Yg'), is given by 

2 

Sty - 2^Sx - 2^(3F C ^ + 2Y c 2 g' ! !)6y , (53) 
a c y c 

where A c = [g c + 5Y c g' c + 2Y?g'l] ~\ 

Let us study linear perturbations about a scalar-field 
dominant solution characterized by Eq. (|35|l and a scaling 
solution characterized by Eq. (|36fl . We recall that these 
critical points satisfy the equations (|30|) and 131|) . Then 
from Eqs. (|28J) and (|29(l we get the following perturbation 
equations: 



d / Sx 



= M 



5x 
Sy 



dN \ Sy 

where the elements of the matrix A4 arc 



(54) 



an 



a-12 



-3 + ^-(2Q + X)x c + 3x 2 c (g c + Y c g' c ) 
-3w m A c x c , 



-3g' c x c Y? 



(55) 



3x c 



V6(Q + X)Y C 



+3w m x c Y c (A c 1 - g c - 2Y c g' c ) 

+v - 6Ac (Q + m, 



Wc 



(56) 



Vc 

021 = y 



-VQX + 6x c {(1 - w m )g c + (1 - 5w m )Y c g' c } 



a-22 



-12w m x c Y c z g' c 
3g' c x 2 c Y c (3w m - 1) + Qw m x 2 c Y 2 9c ■ 



(57) 
(58) 



We checked that these agree what was obtained in 
Ref. [3lJ for three dark energy models. 
The eigenvalues of the matrix M. are 



an + a22 



1±4 1 



4(a 11 a 2 2 - 012021) 
(an + a 2 2) 2 



• (59) 



The fixed point is stable against perturbations when both 
/i + and /i_ are negative or when /i± have negative real 
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parts. Meanwhile it is unstable when either /i + or /i_ is 
positive. 

In what follows we study the case of a non-relativistic 
barotropic fluid (w m = 0). We shall consider the stability 
of fixed points for (i) the scalar-field dominant solution 
and (ii) the scaling solution, separately. 



A. Scalar- field dominant solutions 

When ft,/, = 1 we have (g c + Y c g' c )x c = from 
Eq. (3TJ). Then we find a 2 i = from Eq. This 
means that the eigenvalues of the matrix M. are 



jU+ = an 



/Lt_ = a 22 = -3 + ~^-^ x c 



V6(Q + X)x c , 
V6, 



(60) 
(61) 



Making use of Eq. (|37[) . the second eigenvalue is expressed 
as /j_ = — (3/2)(l — u^). Hence one has \i- < for 
W0 < 1. The first eigenvalue [i+ is negative for 



V6 



2(Q + A) 



(62) 



Since we are considering the case of positive values of Q 
and A, \x- becomes automatically negative when the con- 
dition l|62|) is satisfied. Hence the scalar-field dominant 
fixed point (Q^, = 1) is stable for x c < y/E/2(Q + A), 
whereas it is unstable for x c > v / 6/2(<5 + A). 

In the case of a non-phantom scalar field characterized 
by Px > 0, the stability condition l|62|) can be written as 



Px > 



A(Q + A) 



(63) 



where we used pxx c = A/V6. Meanwhile for a phantom 
field this stability condition corresponds to px < X(Q + 
A)/3. Since px < this is automatically satisfied for 
positive values of Q and A and also for Q = with any 
A. Hence the scalar-field dominant solution is stable for 
phantom fields. This was found in Rcf. |3l| for several 
models of dark energy, but we showed that this property 
persists for any form of g{Y). 



B. Scaling solutions 

The scaling solution satisfies Eq. (|44|l . When w m = 
we have the relation g c x 2 = —Q/(Q + A) and also 
^0 = -Q/(Q + A) + 2Y c g' c x 2 c . Hence we find that the 
components of the matrix A4 arc 



ai2 



a-21 



-an 



0,22 



-3Y c g' c x 2 c . 



-a 2 2 , 



(64) 



From Eq. (|59() we get the eigenvalues of the matrix M. , 



as 



where 



i± yr - ^ 



ii = 

6 = 



3(2Q + A) 
4(Q + A) ' 

(1 -» )(Q + A) 3 [» (g + A) +Q] 
3(2Q + A) 2 



(65) 



(66) 



A c . (67) 



This shows that £i < for positive Q and A. When 
Q = 0, £i is automatically negative (£i = —3/4). /i + is 
negative or has a negative real part. Meanwhile the sign 
of /!_ depends upon that of £ 2 - When £ 2 > the scaling 
solution is stable, whereas it is unstable for £ 2 < 0. 

In order to get a viable scaling solution we require the 
condition fi^ < 1. Hence the stability of the scaling 
solution is dependent on the sign of A c , which can be 
expressed as 



A c = (p Xc +2X c px e x a ) 



This quantity is related to the speed of sound: 



Px 



Px Px + 2Xp X x 



(68) 



(69) 



which appears as a coefficient of the k 2 /a 2 term in pertur- 
bation equations [4lLl42j (fc is a comoving wavenumber). 

While the classical fluctuations may be regarded to be 
stable when c 2 > 0, it was shown in Ref. that the 
stability of quantum fluctuations requires the conditions 
px > and px + 2Xpxx > 0. Hence if we impose 
the quantum stability, A c is positive. Then the scaling 
solution with Sl^ < 1 is a stable attractor. In the case 
of an ordinary scalar field with p = X — V(<f) one has 
A c = 1, which shows that the scaling solution is stable. 
Meanwhile for an ordinary phantom field with p = —X — 
V ((/)), it is unstable since A c = — 1. 

From Eq. I|45|l with w m = we find that the condition, 
< 1, corresponds to 



Px < 



A(Q + A) 



(70) 



When this is satisfied, the stability condition l|63|) for the 
scalar-field dominant fixed point is violated. Hence when 
we have a stable scaling solution with < 1 and A c > 0, 
the scalar-field dominant fixed point with px > is al- 
ways unstable. Hence the system chooses either a scaling 
solution with £7^ < 1 or a scalar-field dominant solution 
with f2<£ = 1 and px > as an attractor. In Ref. j^H this 
property was found by specifying the functional form of 
g{Y), but we have shown that this holds for all scalar- 
field models which possess scaling solutions. 

We note that the scaling solution can be an attractor 
even for a phantom (px < 0) provided that the condi- 
tion px + 2Xpxx > is satisfied. In the case of the 
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dilatonic ghost condensate model (p = —X + ce x ^X 2 ), 
for example, this is realized for 1/3 < cY < 1/2. The 
scalar-field dominant solution cY\ = 1/2 + A 2 /_(A)/16 in 
Eq. (|42|l also exists in the phantom region characterized 
by 1/3 < cYi < 1/2. In the subsection A we showed that 
this case is also a stable attractor. We recall that x\ in 
Eq. (|42|) is negative, whereas x in Eq. I|44|l is positive. We 
have numerically found that the solutions approach the 
scaling solution for the initial condition xi > whereas 
they approach the phantom scalar-field dominant point 
when Xi < 0. Note that another scalar-field dominant 
solution cY 2 = 1/2 + A 2 /+(A)/16 (> 1/2) in Eq. g2J> is 
unstable when the scaling solution is stable. We caution 
that the phantom field {px < 0) with px + 2Xpxx > 
corresponds to a negative sound of speed, thus suffering 
from the instability of perturbations even at the classical 
level 0. See Refs. ElE! for related works. 



V. MULTI-FIELD DARK ENERGY 

We can generalize the analysis to the case in which 
many fields are present with a same Lagrangian density 
for the each field. In the case of an ordinary scalar field 
with an exponential potential several fields can cooper- 
ate to drive an accelerated expansion even if none is able 
to do individually [l^, EH, 01 • We will show that this 
generally happens for the Lagrangian density which pos- 
sesses cosmological scaling solutions. We shall study the 
case without a coupling between dark energy and dark 
matter (Q = 0) throughout this section. This is because 
an accelerated expansion is realized even for Q = in the 
presence of multiple scalar fields. 

Let us consider n scalar fields (<f>i, <f>2, • • • , </>n) with the 
Lagrangian density: 



J2x i9 (x ie x ^), 



(71) 



where Xi = —g^d^idv^ijl and g is an arbitrary func- 
tion. In the presence of a barotropic fluid with an equa- 
tion of state w m = Pm/Pm, the constraint equation is 
given by Cl m + Sl^ = 1 with a fractional density 



(72) 



i=l 



We also obtain the equation for the Hubble rate: 



1 dH 
HdN 



- 1) 



i=l 



.(73) 



Each scalar field, 
(Pjl with Q = 0. 



i>i, satisfies the equation of motion 
Defining dimcnsionless quantities 



l /V&H and y, l = e: 



Ai0i/2 



, we get the following 



differential equations: 



dxj 
dN 



dN 



3xj 
2 



1 + Y j g{Y i )xf - io m (O - 1) — XiXt 



i=l 



+ ^ [A^ - V6{g(Yi) + Y i g , (Y i )}x^ , (74) 



3yi 

2 



^/6 

X? - Wm(&<t> - 1) ^-\Xi 



(75) 



Then the fixed points we are interested in (j/j ^ 0) satisfy 
yf&igiY^+Yig'iYi)] 



g{Yi) + 2Y i g'(Y i ) 
/fi 



(76) 



(77) 



Since the equation of state of the each scalar field is 
given by = g(Y t )/[g(Y t ) + 2Y i g'(Y i )}, we find 



g{Yi) + 2Y l g'{Y i ) 2 
Then from Eq. (|76[) we get 

Wfc = -1 + —\Xi 



(78) 



(79) 



Equation (|77|l indicates that the quantities A^ are in- 
dependent of i = 1, 2, ■ • ■ ,n. Hence we set 



X\Xi — ■ ■ ■ — Xixi — • • • — X n x n = Xx . 



(80) 



Then from Eq. I|79|) . are same for all scalar fields, 
i.e., 



w <t>i 



(81) 



From Eq. Q76[l the quantities Yi are determined by XiXi, 
which means that Yi are also same independent of i: 



Y\ = ■ ■ ■ = Yi = 
Taking note the relation 



Y„ = Y . 



1=1 



1=1 



we find that Eq. l(77|) yields 

(Q - l)(w0 - w m ) = . 



(82) 



(83) 



(84) 



Hence we have two fixed points: (i) scalar-field domi- 
nant solution with = 1 and (ii) scaling solution with 
= w m . In what follows we shall discuss these cases 
separately. 
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Scalar- field dominant solutions (fi^ = 1) and 
assisted dark energy 



From Eqs. JBDJ and (JH2J, Eq. JEJ is written as 
VQ[g(Y) + Yg'(Y)} 



Ax 

Using Eq. (JHJJ with fl 
by 



(85) 



<?(y) + 2Yy(y) • 

1, wc find that Eq. O is given 



Ax = 



Hence if wc choose 



n 1 

i + 9 (y)x 2 A 2 £^ 



1 - 1 

A2 = E I? ' 



(86) 



(87) 



Equation (fH^|) yields 



Equations H85|) and (|88|l show that the system effec- 
tively reduces to that of the single field with A given 
by Eq. l|S7|) . From Eq. ijg5|> and = 1 the quantity 
px = g{y) + Yg'(Y) satisfies the relation s/6xpx = A. 
Hence from Eq. (J75J) with Eqs. (JSOJ and (JHU, we get the 
equation of state: 



W0 



-1 



A 2 
3px 



(89) 



which is the same expression as Eq. (|39|) . 

Let us consider the case of a non-phantom scalar field 
(px > 0). Equation 18 7|) shows that the presence of 
multiple scalar fields leads to the decrease of A 2 relative to 
the single-field case. The quantity Y is known by Eq. H85|) 
in terms of Ax once g(Y) is specified. Then px(Y) = 
g(Y) + Yg'(Y) is determined by the quantity Ax which is 
unaffected by adding scalar fields. This means that the 
effect of multiple scalar fields to the equation of state only 
appears for A 2 in Eq. (|89|) . This effect works to reduce the 
equation of state toward w$ = — 1 as A decreases. Hence 
even if inflation does not occur in the single-field case 
because of large A, the presence of many scalar fields can 
lead to an accelerated expansion by reducing A. Thus we 
have shown that assisted acceleration occurs for all dark 
energy models which possess scaling solutions. 

In the case of a phantom field (px < 0) the presence of 
many scalar fields leads to the increase of w^ toward — 1, 
which is different from the assisted acceleration. In any 
case the equation of state approaches that of cosmological 
constant (w$ = —1) as we add more scalar fields. 

Taking note the relation g{Yi)xl = w^Clfc and 
g(Y)x 2 = u>0 together with l|81|l and (|82[1 . the fractional 
energy density of each scalar field is 



A 2 
A 2 



(90) 



In the two-field case, for example, one has = 
A^/(A 2 + A|) andft 02 = A 2 /(A 2 + A^). The energy density 
is distributed so that the condition l|8U[) is satisfied. 



B. Scaling solutions 

The scaling solution satisfies w ! f >i = w$ = w rn (i = 
1, 2, ■ ■ • , n). Hence we do not have an accelerated expan- 
sion for a non-relativistic dark matter (w m = 0). From 
Eq. (ZZJ and (JHHJ we obtain 



30 i 



V&il + Wrn) 



2X, 



(91) 



This is similar to the single field case, see Eq. (1441 with 
Q = 0. From Eq. (|7T)|) one gets = y/6xipx / '\- Then 
by using Eq. (J5TJ we find Cl^ = 3(1 + u> m )px/A 2 . This 
gives the total energy fraction of the scalar field: 



3(1 



w, 



x 



A 2 



(92) 



where A is defined in Eq. 18 7|) . In the case of an ordinary 
scalar field withpx = 1, this result agrees with what was 
obtained in Ref. [34j . 

In order for the scaling solution to be physically mean- 
ingful, wc require the condition < 1, i.e., 



Px < 



3(1 + W rn ) 



(93) 



In the case of a non-relativistic dark matter this is equiv- 
alent to the condition (|70|l with Q — 0. Since the multi- 
field system is now reduced to an effectively single-field 
system, the stability analysis of critical points discussed 
in Sec. lIVI persists in the presence case as well. Thus the 
scaling solution is stable if a scalar-field dominant solu- 
tion with px > is unstable, and vice versa. This is 
actually what was found in Ref. |34j for an ordinary field 
with an exponential potential. 

As discussed in Ref. [34| we can consider an interest- 
ing situation in which the solution is in the scaling regime 
characterized by w^ >i = w m during radiation and matter 
dominant eras and then enters the regime of an acceler- 
ated expansion as more and more fields join the assisted 
scalar-field dominant attractor. Although this requires 
a fine-tuning, it is interesting that the presence of many 
scalar fields leads to the possibility of an accelerated ex- 
pansion. 



VI. CONCLUSIONS 

In this paper we have provided a general method to de- 
rive critical points relevant to dark energy for the models 
which possess scaling solutions. The existence of scaling 
solutions restricts the Lagrangian density to be in the 
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form p = Xg(Xe x ' t ') where g is an arbitrary function. 
This includes a wide variety of dark energy models such 
as quintessence, phantoms, tachyons and dilatonic ghost 
condensates. Since A is constant, we obtain autonomous 
equations (|28|l and (|29|l for two variables x and y defined 
in Eq. 1|23|) . For the models in which scaling solutions 
do not exist, one has another differential equation for A 
since A is a dynamically varying quantity. Even in this 
case we can apply the results obtained in this paper by 
considering "instantaneous" critical points (x(N),y(N)) 

M 

The fixed points relevant to dark energy is (i) the 
scalar-field dominant solution (f^ = 1) and (ii) the scal- 
ing solution with p^j p m = constant ^ 0. In the former 
case the scalar-field equation of state is generally given 
by Eq. I|37|) . The fixed points can be derived by using 
Eq. (|40() when we specify the function g(Y). Depending 
upon the models of dark energy the fixed points exist 
in the region given by px > or px < 0. In some of 
the models like dilatonic ghost condensate there are two 
fixed points both in the regions p x > and px < 0. 
In the case of scaling solutions the variable x and the 
effective equation of state w e s are uniquely determined 
as Eqs. (I44|) and (|18|l independently of the models. The 
variables y and Y = x 2 /y 2 are known by Eq. (|46|l once 
we specify the form of g(Y). 

The stability of fixed points was discussed in Sec. II VI in 
the presence of a non-relativistic barotropic fluid without 
specifying any form of g(Y). The stability condition for 
the scalar-field dominant solution is given by Eq. (|03l) for 
a non- phantom field (px > 0). When this condition is 
satisfied, the scaling solution exists in the region which is 
not physically meaningful (Q^ > 1). When the quantity 
A c defined in Eq. (|68[l is positive, the scaling solution 
is a stable attractor for 0^ < 1. The positivity of A c 
is ensured if we impose the stability of quantum fluc- 
tuations. Hence when the stable scaling solution exists 



in the region £7^ < 1, the scalar-field dominant solution 
with px > is unstable. 

The above discussion shows that the solution chooses 
cither the scalar-field dominant solution with px > or 
the scaling solution with 0,$ < 1. When there exists a 
phantom scalar-field dominant fixed point (px < 0), we 
found that this solution is classically stable. If both the 
scaling solution and this phantom scalar-field dominant 
solution arc stable critical points, the solutions choose 
either of them depending upon initial conditions. 

Finally we generalized our analysis to the case in which 
more than one field is present with the same form of the 
Lagrangian density. We showed that the system effec- 
tively reduces to that of the single-field with A modified 
as Eq. l|87(l for the scalar-field dominant solution. Since 
the presence of many scalar fields leads to the decrease 
of A toward 0, the equation of state approaches that of 
cosmological constant as we add more fields. Thus we 
have shown that for a non-phantom fluid the assisted ac- 
celeration occurs for all scalar-field dark energy models 
which have scaling solutions. 

It is really remarkable that the derivation of critical 
points, their stability and assisted acceleration can be 
analyzed in a unified way without specifying any form of 
g(Y). We hope that the results obtained in this paper 
will be useful to constrain scalar-field dark energy models 
from future high-precision observations. It is certainly 
of interest to place strong observational constraints not 
only from SN la but from CMB and LSS by using the 
equations of matter perturbations derived in Refs. |42(. 



ACKNOWLEDGMENTS 

This work is supported by JSPS (Grant No. 30318802). 



[1] A. G. Riess et al, Astron. J. 116, 1009 (1998); Astron. 
J. 117, 707 (1999); S. Perlmutter et al, Astrophys. J. 
517, 565 (1999); J. L. Tonry et al, Astrophys. J. 594, 1 
(2003); R. A. Knop et al, Astrophys. J. 598, 102 (2003); 
A. G. Riess et al, Astrophys. J. 607, 665 (2004). 

[2] D. N. Spergel et al, Astrophys. J. Suppl. 148, 175 (2003). 

[3] M. Tegmark et al, Phys. Rev. D 69, 103501 (2004); 
U. Seljak et al, Phys. Rev. D 71, 103515 (2005). 

[4] D. J. Eisenstein et al, Astrophys. J. 633, 560 (2005); 
C. Blake, D. Parkinson, B. Bassett, K. Glazebrook, 
M. Kunz and R. C. Nichol, Mon. Not. Roy. Astron. Soc. 
365, 255 (2006). 

[5] T. D. Saini, S. Raychaudhury, V. Sahni and 
A. A. Starobinsky, Phys. Rev. Lett. 85, 1162 (2000); 
U. Alam, V. Sahni, T. D. Saini and A. A. Starobinsky, 
Mon. Not. Roy. Astron. Soc. 354, 275 (2004); A. Mel- 
chiorri, L. Mersini-Houghton, C. J. Odman and M. Trod- 
den, Phys. Rev. D 68, 043509 (2003); J. Weller and 
A. M. Lewis, Mon. Not. Roy. Astron. Soc. 346, 987 



(2003); T. Padmanabhan and T. R. Choudhury, Mon. 
Not. Roy. Astron. Soc. 344, 823 (2003); T. R. Choud- 
hury and T. Padmanabhan, Astron. Astrophys. 429, 807 
(2005); E. V. Linder and A. Jenkins, Mon. Not. Roy. As- 
tron. Soc. 346, 573 (2003); Y. Wang and M. Tegmark, 
Phys. Rev. Lett. 92, 241302 (2004); Y. Wang and 
P. Mukherjee, Astrophys. J. 606, 654 (2004); S. Nesseris 
and L. Perivolaropoulos, Phys. Rev. D 70, 043531 (2004); 
P. S. Corasaniti, M. Kunz, D. Parkinson, E. J. Copeland 
and B. A. Bassett, Phys. Rev. D 70, 083006 (2004); 
B. A. Bassett, P. S. Corasaniti and M. Kunz, Astro- 
phys. J. 617, LI (2004); B. Feng, X. L. Wang and 
X. M. Zhang, Phys. Lett. B 607, 35 (2005); D. Huterer 
and A. Cooray, Phys. Rev. D 71, 023506 (2005); 
R. Lazkoz, S. Nesseris and L. Perivolaropoulos, JCAP 
0511, 010 (2005); M. Sahlen, A. R. Liddle and D. Parkin- 
son, Phys. Rev. D 72, 083511 (2005); Y. G. Gong, Class. 
Quant. Grav. 22, 2121 (2005). 
[6] V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. D 9, 



10 



373 (20 00); V. S ahni, Lec t. Notes Phys. 653, 141 (2004) 

larXiv:astro-p h/0403324| . 
[7] T. Padmanabhan, Phys. Rept. 380. 235 (2003). 
[8] P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 

(2003). 

[9] E. J. Copeland, M. Sami and S. Tsujikawa, 
arXiv:hep-th/0603057 
[10] A. Y. Kamenshchik, U. Moschella and V. Pasquier, Phys. 
Lett. B 511, 265 (2001); M. C. Bento, O. Bertolami and 

A. A. Sen, Phys. Rev. D 66, 043507 (2002). 

[11] G. R. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B 
485, 208 (2000); V. Sahni and Y. Shtanov, JCAP 0311, 
014 (2003). 

[12] S. C apozziello , S. Carloni and A. Troisi, 
|arXiv:astro-ph /0303041 S. Capozziello, V. F. Car- 
done, S. Carloni and A. Troisi, Int. J. Mod. Phys. D 12, 
1969 (2003); S. M. Carroll, V. Duvvuri, M. Trodden and 
M. S. Turner, Phys. Rev. D 70, 043528 (2004); S. Nojiri 
and S. D. Odintsov, Phys. Rev. D 68, 123512 (2003); 
S. Capozziello, V. F. Cardone, E. Elizalde, S. Nojiri and 
S. D. Odintsov, Phys. Rev. D 73, 043512 (2006). 

[13] Y. Fujii, Phys. Rev. D 26, 2580 (1982); L. H. Ford, 
Phys. Rev. D 35, 2339 (1987); C. Wetterich, Nucl. Phys 

B. 302, 668 (1988); B. Ratra and J. Peebles, Phys. 
Rev D 37, 321 (1988); Y. Fujii and T. Nishioka, Phys. 
Rev. D 42, 361 (1990); E. J. Copeland, A. R. Lid- 
die, and D. Wands, Ann. N. Y. Acad. Sci. 688, 647 
(1993); R. R. Caldwell, R. Dave and P. J. Steinhardt, 
Phys. Rev. Lett. 80, 1582 (1998); I. Zlatev, L. M. Wang 
and P. J. Steinhardt, Phys. Rev. Lett. 82, 896 (1999); 
P. J. Steinhardt, L. M. Wang and I. Zlatev, Phys. Rev. 
D 59, 123504 (1999). 

[14] T. Chiba, T. Okabe and M. Yamaguchi, Phys. Rev. D 62, 
023511 (2000); C. Armendariz-Picon, V. Mukhanov, and 
P. J. Steinhardt, Phys. Rev. Lett. 85, 4438 (2000); C. 
Armendariz-Picon, V. Mukhanov, and P. J. Steinhardt, 
Phys. Rev. D 63, 103510 (2001). 

[15] G. W. Gibbons, Phys. Lett. B 537, 1 (2002); T. Pad- 
manabhan, Phys. Rev. D 66, 021301 (2002); J. S. Bagla, 
H. K. Jassal and T. Padmanabhan, Phys. Rev. D 67, 
063504 (2003). 

[16] R. R. Caldwell, Phys. Lett. B 545, 23-29 (2002); 
R. R. Caldwell, M. Kamionkowski and N. N. Wein- 
berg, Phys. Rev. Lett. 91, 071301 (2003); S. M. Carroll, 
M. Hoffman and M. Trodden, Phys. Rev. D 68, 023509 
(2003); P. Singh, M. Sami and N. Dadhich, Phys. Rev. 
D 68 023522 (2003). 

[17] N. Arkani-Hamed, P. Creminelli, S. Mukohyama and M. 
Zaldarriaga, JCAP 0404, 001 (2004). 

[18] F. Piazza and S. Tsujikawa, JCAP 0407, 004 (2004). 

[19] P. G. Ferreira and M. Joyce, Phys. Rev. Lett. 79, 4740 
(1997); Phys. Rev. D 58, 023503 (1998). 

[20] E. J. Copeland, A. R. Liddle and D. Wands, Phys. Rev. 
D 57, 4686 (1998). 

[21] A. P. Billyard, A. A. Coley and R. J. van den Hoogen, 
Phys. Rev. D 58, 123501 (1998); A. R. Liddle and 
R. J. Scherrer, Phys. Rev. D 59, 023509 (1999); 
J. P. Uzan, Phys. Rev. D 59, 123510 (1999); R. J. van 
den Hoogen, A. A. Coley and D. Wands, Class. Quant. 
Grav. 16, 1843 (1999); A. Nunes and J. P. Mimoso, Phys. 
Lett. B 488, 423 (2000); C. Rubano and J. D. Barrow, 
Phys. Rev. D 64, 127301 (2001); Z. K. Guo, Y. S. Piao, 
R. G. Cai and Y. Z. Zhang, Phys. Lett. B 576, 12 (2003); 
Z. K. Guo, Y. S. Piao and Y. Z. Zhang, Phys. Lett. B 568, 



1 (2003); S. Mizuno, S. J. Lee and E. J. Copeland, Phys. 

Rev. D 70, 043525 (2004); E. J. Copeland, S. J. Lee, 

J. E. Lidsey and S. Mizuno, Phys. Rev. D 71, 023526 

(2005); M. Sami, N. Savchenko and A. Toporensky, Phys. 

Rev. D 70, 123526 (2004); V. Pettorino, C. Baccigalupi 

and F. Perrotta, JCAP 0512, 003 (2005). 
[22] A. de la Macorra and G. Piccinelli, Phys. Rev. D 61, 

123503 (2000); S. C. C. Ng, N. J. Nunes and F. Rosati, 

Phys. Rev. D 64, 083510 (2001). 
[23] T. Barreiro, E. J. Copeland and N. J. Nunes, Phys. Rev. 

D 61, 127301 (2000). 
[24] V. Sahni and L. M. Wang, Phys. Rev. D 62, 103517 

(2000). 

[25] L. Amendola, Phys. Rev. D 62, 043511 (2000); L. Amen- 

dola, Phys. Rev. D 60, 043501 (1999). 
[26] L. Amendola and D. Tocchini-Valentini, Phys. Rev. D 64, 

043509 (2001); D. Tocchini-Valentini and L. Amendola, 

Phys. Rev. D 65, 063508 (2002). 
[27] L. R. W. Abramo and F. Finelli, Phys. Lett. B 575, 165 

(2003). 

[28] J. M. Aguirregabiria and R. Lazkoz, Phys. Rev. D 69, 
123502 (2004). 

[29] E. J. Copeland, M. R. Garousi, M. Sami and S. Tsu- 
jikawa, Phys. Rev. D 71, 043003 (2005). 

[30] S. Tsujikawa and M. Sami, Phys. Lett. B 603, 113 (2004). 

[31] B. Gumjudpai, T. Naskar, M. Sami and S. Tsujikawa, 
JCAP 0506, 007 (2005). 

[32] A. R. Liddle, A. Mazumdar and F. E. Schunck, Phys. 
Rev. D 58, 061301 (1998). 

[33] A. A. Coley and R. J. van den Hoogen, Phys. Rev. D 62, 
023517 (2000). 

[34] S. A. Kim, A. R. Liddle and S. Tsujikawa, Phys. Rev. D 
72, 043506 (2005). 

[35] M. Doran and J. Jaeckel, Phys. Rev. D 66, 043519 (2002); 
L. P. Chimento, A. S. Jakubi, D. Pavon and W. Zimdahl, 
Phys. Rev. D 67, 083513 (2003); M. Sami and T. Pad- 
manabhan, Phys. Rev. D 67, 083509 (2003); L. Amen- 
dola and C. Quercellini, Phys. Rev. D 68, 023514 (2003); 
M. B. Hoffman, arXiv:astro-ph/0307350 A. V. Maccio 
et al, Phys. Re v. D 69, 1235 16 (2004); G. Huey and 
B. D. Wandelt, |arXiv:astro-ph/0407196| L. Amendola, 
M. Gasperini and F. Piazza, ' |arXlv!astro-p h/0407573 
S. Nojiri, S. D. Odintsov and S. Tsujikawa, Phys. Rev. 
D 71, 063004 (2005); S. Nojiri and S. D. Odintsov, 
Phys. Rev. D 72, 023003 (2005); G. Olivares, F. Atrio- 
Barandela and D. Pavon, Phys. Rev. D 71, 063523 
(2005); G. Calcagni, S. Tsujikawa and M. Sami, Class. 
Quant. Grav. 22, 3977 (2005); T. Koivisto, Phys. Rev. 
D 72, 043516 (2005); B. Wa ng, C. Y. Lin and E. Ab- 
dalla, |arXiv:hep-th/ 0509107 S. Lee, G. C. Liu and 
K. W. Ng,|arXiv:astro~ /0601333 B. Hu and Y. Ling, 
|arXiv:hep-th/0601093| 

[36] J. P. Uzan, Phys. Rev. D 59, 123510 (1999); T. Chiba, 
Phys. Rev. D 60, 083508 (1999); N. Bartolo and 
M. Pietroni, Phys. Rev. D 61 023518 (2000); F. Per- 
rotta, C. Baccigalupi and S. Matarrese, Phys. Rev. D 
61, 023507 (2000); B. Boisseau, G. Esposito-Farese, 
D. Polarski and A. A. Starobinsky, Phys. Rev. Lett. 85, 
2236 (2000); G. Esposito-Farese and D. Polarski, Phys. 
Rev. D 63, 063504 (2001); E. Elizalde, S. Nojiri and 
S. D. Odintsov, Phys. Rev. D 70, 043539 (2004). 

[37] P. Gu, X. Wang and X. Zhang, Phys. Rev. D 68, 
087301 (2003); R. Fardon, A. E. Nelson and N. Weiner, 
JCAP 0410, 005 (2004); X.-J. Bi, P. Gu, X. Wang and 



11 



X. Zhang, Phys. Rev. D 69 113007 (2004); P. Q. Hung 
and H. Pas, Mod. Phys. Lett. A 20, 1209 (2005); 
D. B. Kaplhan, A. E. Nelson and N. Weiner, Phys. 
Rev. Lett. 93, 091801 (2004); R. D. Peccei, Phys. 
Rev. D 71 02352 7 (2005); E. I Guendelman and 
A. B. Kaganovich, |hep-th/04lTT88] X. J. Bi, B. Feng, 
H. Li and X. m. Zhang, Phys. Rev. D 72, 123523 (2005); 
A. W. Brookfield, C. van de Bru ck, D. F. Mota and 
D. Tocchini-Valentini, |astro-ph /0503349 A. W. Brook- 
field, C. van de Bruck, D. F. Mota and D. Tocchini- 
Valentini, |arXiv:astro-ph/ 0512367 | H. Li, B. Feng 
J. Q. Xia and X. Zhang, |arXiv:astro^p h/0509272 
P. H. Gu, X. J. Bi, B. Feng, B. L. Young and X. Zhang, 
|arXiv:hep-ph/0512076| 



[38] S. Tsujikawa, Phys. Rev. D 72, 083512 (2005). 

[39] D. Comelli, M. Pietroni and A. Riotto, Phys. Lett. B 

571, 115 (2003). 
[40] L. Amendola, M. Quartin, I. Waga and S. Tsujikawa, 

work in progress. 
[41] J. Garriga and V. F. Mukhanov, Phys. Lett. B 458, 219 

(1999); J. c. Hwang and H. Noh, Phys. Rev. D 66, 084009 

(2002). 

[42] L. Amendola, Phys. Rev. Lett. 93, 181102 (2004); 

L. Amendola, S. Tsujikawa and M. Sami, Phys. Lett. 

B 632, 155 (2006). 
[43] A. Vikman, Phys. Rev. D 71, 023515 (2005); A. Anisimov 

and A. Vikman, JCAP 0504, 009 (2005). 



